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ESTIMATING SURFACE TEMPERATURE IN' FORCED CONVECTION 
NUCLEATE BOILING - A SIMPLIFIED METHOD 
by R. C. Hendricks and S. S. Papell 

NASA Lewis Research Center 
Cleveland, Ohio 

ABSTRACT 

A simplified expression to estimate surface temperatures in forced 
convection boiling was developed using a liquid nitrogen data base. Using 
the principal of corresponding states and the Kutateladze relation for max- 
imum pool boiling heat flux, the expression was normalized for use with 
other fluids. The expression was applied also to neon and water. For 
the neon data base, the agreement was acceptable with the exclusion of 
one set suspected to be in the transition boiling regime. For the water 
data base at reduced pressure greater than 0.05 the agreement is generally 
good. At lower reduced pressures, the water data scatter and the calculated 
temperature becomes a function of flow rate. 


The estimation of wall temperature in forced convection boiling 
usually involves lengthy calculations. However, a simplified empirical 
technique was developed using the liquid nitrogen data base of ref. 1, 
the principle of corresponding states and the Kutateladze expression for 
maximum nucleate heat flux in pool boiling ref. 2. 

The empirical correlation in reduced form is 
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where 




p/p c 
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a; = Pitzer's accentric factor 
q r = ^/ q k 

The critical pool heat flux 

q*=°. 16h fg (gcrp|p fg )1/4 (2) 

was deduced from the Kutateladze expression ref. 2, and occurs for 
most fluids near P r = 0. 3. Values of for some fluids and a method 
to estimate when data are not available are given in Table I. 

The data set for forced convection boiling nitrogen ref, 1 was ex- 
amined for developed boiling cases. While the wall temperature peaks 
at incipience and decreases downstream of this point, the difference is 
not large and the surface temperature 10 to 20 diameters downstream 
was chosen as the forced convection boiling surface temperature, T 

w 

Knowing the pressure at this point defines the saturation conditions and 
q£, and the critical constants follow from Table I. Figure 1 displays 
the calculated temperature difference (T - T ) as a function of the ex- 
perimental value for nitrogen. The agreement is generally good. 

With an equation established, it is necessary to determine if it can 
be applied to other nonpolar fluids. Using the forced convection boiling 


3 


neon data of ref. 1, and following the previously defined procedure, the 
calculated and experimental temperature differences are displayed in 
fig. 2. Again the agreement is reasonably good, with the exception of 
one data set at 0. 7 MPa. The level appears to be about IK higher than 
calculated. Subsequent examination of the rteon data indicates that tran- 
sition boiling can take place (followed by film boiling - see point on fig- 
ure 2) and it is suspected that these data may be near the transition boiling 
regime. 

A rather difficult test of the equation would be to apply it to a polar 

molecule, e. g. , water. Sources of forced convection boiling water data 

such as McDougal, ref. 3; Schaffer, ref. 4; UCLA, ref. 5; Clark and 

Rohsenow, ref. 6; and Kreith and Summerfield, ref. 7; were considered 

by Papell, ref. 8. Selecting a few of these points, the calculated and 

experimental temperature differences are given in figure 3. Here the 

agreement is not as good as for the nonpolar molecules (especially for 

P r < 0. 05) but for the most part the agreement is reasonable considering 

the large critical temperature and pressure of water. Reference 5 points 

out that at low pressure, bubble mechanics have an increased influence on 

the flow rate which, in turn, makes T r somewhat dependent on flow rate. 

For water at pressures less than 0. 1 P this technique may predict 

(T - T ) only within 50 percent but in most cases it is conservative, 
w s 

NOTATION 

g gravitational acceleration (cm/sec ) 

h„ heat of vaporation (j/g) 

S 2 
P pressure (MN/m ) 

P reduced pressure (P/P ) 

X L 

2 

q heat flux (w/ cm ) 
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Kutateladze critical pool heat flux (w/cm ) 

gas constant . -1- 

g K 

temperature, K 
temperature difference, K 

3 

specific volume (cm /g) 
compressibility 

O 

density (g/cm ) 
surface tension 

Pitzers accentric factor (see ref. 9) 

Subscripts 
c critical 

fg liquid to vapor 

g vapor 

r reduced 

s saturated 

w wall 
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TABLE I. - q max FROM THE KUTATELADZE POOL BOILING RELATION 
AND ASSOCIATED PARAMETERS 


Fluid 

T 

c 

K 

P c 2 
MN/m‘ 

p c 

g/ec 

Z c 

Kutateladze 
maximum - heat flux 






T P q 

K 

MN/m 2 

W/cm 2 

Helium 

5.2 

0. 2275 

0.0693 

0.304 

4.06 

0.086 

0.77 

Para-Hydrogen 

32.576 

1.291 

.03143 

.302 

26. 8 

.409 

12.4 

Noon 

44. 4 

2.654 

.483 

.300 

36.3 

.796 

19. 1 

Nitrogen 

126.3 

3.417 

,3105 

.294 

104 

1.04 

34 

Carbon 

132, 01 

3. 499 

.2997 

. 296 

no 

1.07 

35. 4 

Monoxide 








Fluorine 

144.31 

5.215 

.5738 

.288 

123 

1, 89 

6 * 

Argon 

150.7 

4. 865 

. 531 

.292 

124 

1.52 

4-,. 3 

Oxygen 

154. 78 

5.082 

. 4325 

. 202 

128 

1.55 

49. 8 

Methane 

190.77 

4.626 

,162 

. 289 

157 

1.39 

58.5 

Carbon 

304.21 

7.383 

,464 

.277 

257 

2.22 

98.3 

Dioxide 








Water 

647. 2 

22.08 

.32 

. 233 

558 

6.9 

489.3 


q mnx may be approximated as: q max = — (Pg)’' ^ where P c in MN/m‘ anti 
q is W/cm 2 to +25 percent which is equivalent to +0 percent In AT„ 


Figure I. - Accuracy of correlation equation (1) lor predicting surface 
lemperature - LIQUID NITROGEN DATA. 
















